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ABSTRACT 



The ability to predict torsional instability in the 
early stages of design can have important consequences 
on the design of both conventional and high performance 
ships. This thesis develops fast approximate methods of 
torsional buckling analysis for particular application 
in the concept, feasibility, and preliminary stages of 
ship design. 

Two simplified models of stiffeners commonly used in 
ship construction were presented. The first was an 
ideal I section stiffener. The second, a flat bar 
stiffener, which was included in this analysis because 
of its favorable product ional properties and wide 
useage. These models, with initial imperfections, were 
subjected to axial compressive end loading and the 
resultant behavior analyzed. The approach of the 
analysis of this thesis was an application of the 
energy method to determine the critical buckling stress 
instead of the more commonly used equilibrium approach. 
Both beam theory and thin plate theory were used in 
energy-work relat ionships in these derivations. Then 
the first yield load was determined utilizing the 
Perry-Robertson approach. Integrated into all phases of 
this exploration was the concept of the geometric 
imperfect ions of the stiffeners. 
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NOMENCLATURE 



A Nondimensional parameter 

a Length of stiffener between transverse supports 

b Uniform stiffener spacing 

b Plating effective width 

e 

B Angle of rotation of the stiffener about the toe 
* 

B Initial imperfection of the stiffener in terms of a 
rotation angle about the toe 



B Amplitude of the sinusoidal function for B 
o 

* * 

B Amplitude of the sinusoidal function for B 
o 

C Rotational spring constant (moment/length) of the 
supporting plating 

C Longitudinal warping constant of stiffener 
w 

D Flexural rigidity of plating 

D Flexural rigidity of web plating of stiffener 
w 

d Depth of flat bar stiffener; overall depth of Tee 
stiffener 

d Depth of stiffener to midthickness of flange 
c 

d Depth of web 
w 

E Young’s Modulus of the material 



5 



§ 

X 

f 

w 

G 

h 



I 

I 

P 

I 

t 



I 

z 

J 

L 

M(x) 



m 

P 



s 

§> 

CR 



Membrane strain 

Width of stiffener flange 

Shear Modulus of the material 

Height of neutral axis of plate-st if fener 
combination from midplane of plating 

Effective vertical moment of inertia of stiffener 
and associated effective width of plating 

Polar moment of inertia of stiffener about toe 

Vertical moment of inertia of stiffener alone about 
toe 

Moment of inertia of stiffener about web plane 

St. Venant’s torsion constant for stiffener 

Reduced slenderness ratio 

Moment developed in the stiffener by the total 
torque 

Mode number 

Axial end load 

Height of stiffener shear center above toe 

Elastic axial torsional buckling stress in 
stiffener (critical stress) 



£ Elastic 
CL 



axial torsional buckling stress limit 
stiffener (classical stress) 



of 



£ Axial stress in stiffener 
e 



& 



£ Axial stress component 



x 



> Maximum compressive stress developed by 

x MAXT 

the torque 



§ Tensile yield stress of material 



y 



T Amplitude of the sinsoidal torque function 
o 

T Total torque developed in the stiffener structure 
by the loading 

T St. Venant’s torque developed in the stiffener 
SV 



T Warping torque developed in the stiffener 
w 



t Plate thickness 

t Stiffener flange thickness 
f 

t Stiffener web thickness 
w 

u Displacement in the longitudinal (x) direction 
V Total strain energy of the structure 
v Sideways flexure (y) of stiffener 
* 

v Initial imperfection of stiffener in horizontal (y) 
direct ion 

W Total work of the applied force 
w Vertical flexure <z> of stiffener 
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* 

w Initial imperfection of the stiffener in the 
vertical <z) direction 

z Height of stiffener centroid above toe 

■\J Poisson’s Ratio 



Note: x,y or z subscripts indicate partial derivatives 

with respect to those coordinates . These are used 

* * * 
particularly with v, v , w, w , B, B . 
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INTRODUCTION 



The subject of this thesis is to develop a 
simplified analytical method of predicting torsional 
buckling of stiffeners used in ship structures. 
Torsional buckling or instability is the phenomena by 
which a column, at a certain axial load, fails suddenly 
in a combined mode of twist and lateral bending of the 
cross section (reference 7). Figure (0-la) (from 
reference 6) illustrates the failure mode referred to 
here as torsional buckling - in this case for a flat 
bar stiffener. Figure (0-lb) illustrates the failure 
mode of tf*e tripping phenomenon which is not covered in 
this analysis. Figure (0-2) (reference 6) is the 
experimental load deflection curve for the failure of a 
flat bar stiffener under an axial load. Point ft of the 
solid (experimental) curve of figure (0-2) is the 
ultimate strength load. As can be seen from figure (0- 
2), the curve proceeds in a downward fashion - the bar 
looses stiffness and fails from this point. 

The ability to predict torsional instability 
in the early stages of design can have important 
consequences on the design of both conventional and 
high performance ships. It is the intention of this 
thesis to develop fast approximate methods of torsional 
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Figure (0-la) Torsional Buckling 




Figure (0-lb) Tripping 



(0-1) Character i zat ion of Torsional Buckling 
and Tripping in a Flat Bar Stiffener 



e oa 




Figure (0-2) 



Experimental Load Deflection Curve 
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buckling analysis for particular application in the 
concept, feasibility, and preliminary stages of ship 
design. 

To this end, two simplified models of stiffeners 
commonly used in ship construction will be presented. 
The first is an ideal I section stiffener. The second, 
a flat bar stiffener, which is included in this 
analysis because of its favorable product ional 
properties and wide useage. These models, with initial 
imperfect ions, will be subjected to axial compressive 
end loading and the resultant behavior analyzed. The 
approach of the analysis of this thesis will be an 
application of energy methods to^determine the critical 
buckling stress instead of the more commonly used 
equilibrium approach. Both beam theory and thin plate 
theory are used in energy-work relat ionships in these 
derivations. Then the first yield load is determined 
utilizing the Perry-Robertson approach as described in 
references 0,14). Integrated into all phases of this 
exploration will be the concept of the geometric 
imperfections of the stiffeners. 

It is the intention of this thesis, using the above 
mentioned approach, to lay groundwork which would have 
the potential for future expansion to include mode 
interaction - to be able to easily analyze combined 
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loadings etc.. 

Particular emphasis is placed here upon inclusion 
of the geometrical initial deformities of the 
stiffeners which occur during the manufacture and 
through day to day usage of the ship. It is well known 
that these geometric imperfect ions can drastically 
reduce the strength of the structural member. But 
current design standards (references 10,11 etc) have 
not tried to quantify these values. Some design codes 
(example: reference 10) have specified tolerence limits 
of various geometric parameters which must be met 
before the design formulae are considered valid. 
Granted, in the concept, feasibility and preliminary 
stages of ship design where the ship only exists on 
paper, the designer would have no knowledge of the end 
resulting imperfections of the actual ship. Hopefully, 
in the future, with the advent of better production 
control in shipyards and the incursion of the computer 
into ship production facilities, this information will 
be collected and statistically analyzed. 

In addition, once the ship is manufact ured, these 
formulae will provide a quick and easy evaluation of 
the actual critical buckling stress in terms of 
physical parameters. Thus, this analysis has made a 
special effort to define the imperfect ions in easily 
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measured terms. This analysis measures the horizontal 
or vertical movement of the shear center of the 
stiffener at its point of greatest deflection and 
assumes a simple sinusoidal distr i but ion. Knowledge of 
the web depth easily converts the imperfection into a 
rotational angle. Here, the stiffener is considered to 
rotate about its toe or base. In other words, the base 
plating is considered much larger and stiffer. Figure 
(0-3) illustrates this simplified concept for a Tee 
stiffener. 




Figure (0-3) Tee Stiffener 

Currently the method of handling the stiffener 
failure problem is incorporated into the use of factors 
of safety. But, with recent attempts to further 
optimize ship structures - in particular to minimize 



IS 



the amount of weight of a ship that must be devoted to 
its structure - better, yet reliable, methods must be 
devised. Remember, with a ship, less weight in the 
structure equates to more weight for the payload be it 
weapons for a military ship or cargo for a civilian 
ship. 

This thesis concludes with a comparison between the 
formulae derived from the simplified models presented 
here and published formulae, published finite element 
analysis, and exper irnental ly derived results. 
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1. DESCRIPTION OF MODELS 



Torsional buckling or instability is character i ted 
by a twisting of the stiffener about its line of 
attachment to the plating. This deformation pattern 
involves both sideways and vertical flexure (v, w) and 
rotation (B) of the stiffener as shown in figure (1) as 
described in reference (3). 




Figure (1-la) Coordinate System 





Figure (1— lb) Undeformed Figure (1-lc) Deformed 

Figure (1-1) Characteri zat ion of Torsional Buckling 



ia 



is a simpl i f icat ion of 



1. 1 Model Ij_ 

The first model considered 



an ideal I- 


-sect ion. 


It consists 


torsional 


rigidity. 


The web 


of 


only shear 


forces. I 


n effect, 


it 


flanges: 









of an I-beam with zero 
the I-beam transmits 
looks like two unequal 



■ /////’. zjzn 




I-beam Stress Distribution 

Figure (1-2) Model I 

This model has been considered in reference (4). 
However- reference (4) used an equilibrium approach and 
initial imperfections were considered as an 
eccentricity of the axial load. The resulting equations 
of the critical stress were quite complicated. 

Figure (1-3) illustrates the geometric parameters 
of this model. 
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Figure (1-3) - Geometric Torsional Buckling 

Parameters for Model I 



I - Moment of inertia about the web plane (top 
z 

only) 

3 

I = (t f ) / 12 
z f w 

s - Height of shear center above toe (origin) 
flange only) 

s * d 

c 

C - Longitudinal warping constant about shear 

w 

of stiffener alone (top flange 

C = 0 
w 

I - Vertical moment of inertia about toe (top 
t 

only) 

2 2 

I = f t (d + ( 1 / 12) t ) 
t w f c f 



flange 



(top 

center 

only) 

flange 
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Polar moment of inertia about toe 



I - 
P 



1 = I + I 

P t 2 

2 - Height of centroid above toe 

2 = d 



J - St. Venant’s torsion constant (top flange only) 
3 

J = (ft )/3 
w f 



£1 



1.2 Model II.. 

The second model considered is the simple flat bar 
stiffener illustrated in figure (1-4). 




Figure (1—4) Model II - Flat bar stiffener. 



Figure (1-5) illustrates the geometric parameters 
of this model. 
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Figure (1-5) - Geometric Torsional Buckling 

Parameters for Model II 



I - Moment of inertia about the web plane 
2 

(web only) 

3 

I = (d t )/12 
z w w 

s - Height of shear center above toe (origin) 

(web only) 

s = d /2 
w 

C - Longitudinal warping constant about shear 
w 

center of stiffener alone 

C = 0 
w 

I — Vertical moment of inertia about toe (web only) 
t 3 

I « (t d )/3 

t w w 

I - Polar moment of inertia about toe 
P 

I = I + I 

P t 2 
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z - Height of centroid above toe 

z = d /2 
w 

J - St. Venant’s torsion constant 
3 

J = (d t >/3 
w w 

D - Flexural rigidity of the web plate 
w 

3 

Et 

w 

D = 

w 2 

12(1- V ) 
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2. MODEL I - TORSIONAL BUCKLING UNDER AXIAL LOADS 



2. 1 Development of the Strain Energy Eguat jion^ 

Since Model I basically consists of two unequal 
flanges with the top flange rotating about the bottom 
flange, the torsional rigidity is considered only for 
the top flange. Thus the strain energy equation 
consists of three parts (reference 1). The first term 
represents sideways bending. The second term represents 
longitudinal warping. The third term represents the 
torsional rigidity of the top flange alone. The strain 
energy for a length a is: 

r a 2 2 2 

V = 1/2 J (El v + EC B + GJB )dx 

'0 Z XX W XX X 



Also from geometrical properties we have (see figure 
1-3) : 

v — sB 



Thus: 



V 



1/2 



Ct 



a 2 2 2 

(El s B + EC B 

0 Z XX W XX 



/; 



a 2 2 

E(I s + C >8 



2 

+ GJB ) dx 
x 



V 



1/2 



25 



W XX 



2 

+ GJB 3 dx 

x 



The following choice of B is made, 



in order to 



meet 
x=a ) 



certain boundary conditions 
and to keep B simple. 



( ie B 

xx 



0 at x=0. 



B = B sin (mTf'x/a) 
o 

Thus for Model I: 

2 £ 2 £ 

V = 1/4 a B (rnTTVa) C E(I s + C ) (m 7>7a> + GJ 1 

o z w 
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2. 2 Sign Convention 

Before development of the virtual work equation, a 
note should be said about the sign convention used 
here. Compression is held to be a positive stress and 
tension as negative. The following figures illustrate 
the sign convention used here. 




4 







v 



A 



xz 













4 a 



x 



Figure (2-1) Sign convention for inplane axial and 

shear stress components. 
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Figure 




( 2 - 2 ) 



Application of sign 
case of an axi 



convent ion 
al load. 



for 
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2. 3 Development of Virtual. Work Eguat_i on^ 

Model I is not assumed to be perfect but rather to 

have initial imperfections. These initial imperfections 

* * * 

v , w and B and the additional deformations v, w and E 
are illustrated in figure (2-3). 




Figure (2-3) Illustration of initial imperfect ions 
v*, w*, B* and additional deformations v,w, B. 



The virtual work for the case of an inplane axial 

load becomes (reference l)s 

W = J | §> u(y,z) dy dz 

"^A x 

Thus using figure (1-1), it can be shown that the 

strain of the centroidal axis in the bent 

conf igurat ion, & (membrane strain), is (reference 1): 

x 



2 2 

& - -u + l/2(v ) + l/2(w ) 

XX X X 



29 



for a perfect I section. With imperfections, 

subst itute: 

* 

v + v for v 
* 

w +> w for w 



Thus: 

* 2 * £ 

gj = — u + 1 /£ ( v +v ) + l/2(w +w ) 

xx xx xx 

using the inextensibi 1 ity assumptions, set: 

& = 0 
x 

we have: 

*2 * 2 
u = l/2(v +v ) + l/£(w +w ) 

x xx xx 



u = 1 /£ 



L 



a * £ * £ 

C (v ) + (w +w ) 1 dx 

0 x x xx 



where u = u (y, z) . 

From geometrical relat ionsh ips (see figure (£-3)) 



V 


At 


zB 


w 


Ss -yB 


* 




* 


* 


* 


V 


£ 


zB 


w 


* -yB 


V 


At 


zB 


w 


2! -yB 


X 




X 


X 


X 


* 




* 


* 


* 


V 


At 


zB 


w 


2: -yB 


X 




X 


X 


X 



Substituting and integrating we get: 



we 
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/a 2 2 2 * 2 *2 

u (y, z) = 1/2 j (z B +2zBB + z B 



0 



X X 



2 2 2 *2 2 *2 
+ y B + 2y B B +yB >dx 
x xx x 



^a 2 2 * 2 

u (y, z) = 1/2 | (z + y ) (B + B ) dx 

'0 



X X 



Again using the simple choice for B and B : 



B = B sin (m If x/a) 
o 

* * 

B = B sin (m/T x/a) 
o 



Thus:- 



^a 2 2 

u(y, z) = 1/2 j (z +y )C(m7?Va)B cos (m /T x/a) 
y 0 o 



* 2 

+ ((m^/aJB cos(m7Vx/a)3 dx 
o 



ra 2 2 * 2 

u(y, z) = 1/2 J (z +y ) (B +B ) 



0 



O O 



* C <m V' / a) cos (m /^x/a) ] dx 



2 2 * 2 2 
u (y, z) = 1/2 (z +y ) (B +B ) (m/^/a) 

o o 

a 

* Cx/2 + (sin(2m77'x/a)/(4m^'/a)3 

0 
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2 2 * 2 2 
u(y,z) = l/4(z +y ) (B +B ) (rnT^/a) a 

o o 



The work equation now becomes: 

2 2 * 2 
[ 1/4 ( z +y ) (B +E 

ft x 



r r 2 2 * 2 2 

W = JJ £ Cl/4 (z +y ) (B +B ) (m^/a) a3 dy dz 



o o 



For the current situation of end loading, £ = £ 

x e 

* 2 2 rr 2 2 
W = £ C 1/4 (B +B ) <m /77a) a] j J (z +y ) dy d; 
e o o ' ' ft 



But note that: 



// ♦ y 2 ) dy 



dz = I 



Thus: 



* 2 2 

W = 1/4 £ (B +B ) <rn TT'/a) a I 
e o o p 
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2. 4 Determination of the Criticai Buckling Stress.. 

Applying the calculus of variations with respect 

to B to the work and strain energy equations yields 
o 

(reference 5) : 

* 2 2 

6 W = 1/2 £ <B +B ) <m#Va) a I d B 



e o o 



P o 



2 2 2 
cfv = 1/2 a B <m7>Va) CE ( I s + C ) (m D'/a) 



w 



+ GJ1 A B 



S 



Applying the principle of minimum potential energy: 

IT = - £w * 0 

we get : 

S w = 

Solving for £ : 
e 



2 2 
CE(I s + C ) (m D'/a) + GJ] 

z w 



B 



* 

B + B 
o o 



( 2 - 2 ) 

It can be seen by inspection that the lowest 
buckling stress occurs for one wave, m = 1, since m 
must be an integer. Thus: 
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£ = £ 
e CR 

As it can easily be seen, as B becomes 

o 

approaches the limit of: 

2 2 
[E(I s + C ) (m TT/a) + GJ3 
z w 

% 

CL 

I 



large % 

CR 



(2-3) 



P 



Figure (2-4) 
increasing B 

o 



illustrates the behavior of % with 
* CR 

B . 
o 



Further simplification and rearranging of equations 
(2-2) and (2-3) provides: 

. B 

o 

S = £ 

CR CL * 

B + B 
o o 



* 

SB+£B=%B 
CR o CR o CR o 



* 

B 

o 

B 

o 

£ /£ - 1 
CL CR 



(2-4) 
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* 

Figure (£-4) £ vs. (B + B ) 

CR o o 



35 



2.5 Determination of the First Yield Load. 



To determine the first yield load, the total torque 
developed in the torsional buckling process must be 
first evaluated. This torque will consist of two parts. 
The first term consists of the well known St. Venant’s 
Torque (reference 1). The second term represents the 
warping torque (reference 2). 

From reference ( 1 ) : 

T = T + T 

SV W 

where: 

T = GJB 
SV K 

and from reference (2) : 

T = - El d w 

W 2 C XXX 

and from geometry; 

w = d B 
x c x 

resulting in: 

2 

T = - El d B 
W z c xxx 

Thus the total torque is: 
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(2-5) 



2 

T = GJB - El d B 

X Z C XXX 

Figure (2-3) illustrates the above equations. 



Again a simple choice of B is made, and with m = 1: 

3 = B sin(7?-x/a) (2-6) 

o 

Equation (2-5) now becomes: 

T = GJ(7'Va)B cos(// t x/a) 
o 



2 3 

+ El d {77-/ a) B cos(7Tx/a) 

2 C O 



2 2 

T = ()^/a)B cos ( 7^-x/a) CGJ + El d (^/a) 1 



o z c 

Choosing a simple sinusoidal shape for T allows: 

T = T cos {17 x/a) 
o 

Thus: 



2 2 

T = ( V" /a) B CGJ + El d (T/Va) 1 



(2-7) 



z c 



To eliminate B , 



o 

(2-7). 



equation (2-4) is applied to equation 



* 

B 

o 

2 2 

T = CGJ + El d ( 77 ~ / a ) ](7*/ a ) 

o z c 

(£ /% - 1 ) 
CL CR 



(2-3) 
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From sectional considerations, an evaluation of the 



maximum stress can be made. 



rv 

\ 



Figure _(2-5) Sectional illustration. 



The moment developed by the torque in terms of 

deflections is (reference 2): 

M(x) = - El w. 

2 xx 

So the maximum compressive stress developed in the 
flange due to the torque is: 



- El d B f 

2 C XX W 

<S ) 

x MfiXT 

21 



(S > - - Ed B f /2 

x MflXT c xx w 



Application of the above definition for B (equation 
2-4) leaves: 
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Ed f 
c w 



2 

(£ ) = i/T’/a) B sin (77- K /a) 

x MftXT o 



p 



ft choice of x is made such that (S ) is a maximum: 

x MflXT 



Ed f 

c w 2 

(£ ) = ( V'/a) B 

x MftXT o 

2 



(2-9) 



Equations (2-7) and (2-9) are combined to eliminate B . 

o 



Ed f T ( tfVa) 
c w o 

(S ) 

x MftXT 2 2 

2CGJ +EI d ( 7^/a) 1 
2 c 



( 2 - 10 ) 



Using equation (£-3) leads to: 



Ed f <7T/ a) 
c w 



B 



(£ ) 
x MftXT 



CS /S - n 

CL CR 



( 2 - 11 ) 



Note that when B =0 that this equation leads to an 

o 

indeterminate form of 5 since £ = % 

x CL CR 

The total maximum stress can now be described by: 



(§S ) = (£ ) + P/fl 

x TOTfiL x MftXT 



where 
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p/fl = % 



CR 



Thus the yield stress has been shown to be: 

S = (£ ) + % (2-12) 

y x MflXT CR 

Using equation (2-11) in equation (2-12): 

2 * 

Ed f (7 * /a) B 

c w 




2 S /£ - 1 

CL CR 



2 

Solving for £ and dividing by £ to 

CR y 

nondimensional form yields: 



(2-13) 

place into 




(2-14) 



Choosing the smallest root and simplifying the equation 
by introducing a slenderness ratio L: 



% 



L = 




(2-15) 
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and letting 



A 



2 * 
< u /a) Ed f B 
c w o 




(2-16) 



Then s 
£ 

CR 



2 

(1 + 1/L + A) 



N 



(1 + 1/L 



A) 



2 

4/L 



£ 



2 



y 

(2-17) 

A plot of ?» /£ vs L for various values of A is 

CR y 

shown in figure (2-6). Figure (2-6) was compared with 
various design criteria from reference (10) and it was 
found that equation (2-17) was similar to design 
equations of reference (10). However, the reduced 
slenderness ratio, L, of reference (10) is based upon 
simple Euler buckling and the geometric parameter, A, 
is of course different. The column curves of reference 
(10) were quite similar in appearance and range of 
values to that found in figure (2-6). 
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Figure (2-6) Column Curves 
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3. MODEL II - TORSIONAL BUCKLING UNDER AXIAL LOADS 



3. 1 Development of the Strain Energy Eguat ion^ 

Since Model II is the very simple case of a flat 
bar stiffener, thin plate theory can be easily utilized 
in the development of the strain energy equation. Much 
of this section has been patterned after a similar 
analysis presented in reference (3). 

For a web plate of uniform thickness and a length a 
the strain energy is repensented by (reference 5) : 




c. 

- 2(1-N)Cv v - v ]] dx d: 

XX zz xz 

f a ^ 

+ 1/2 I CB dx 
•'o 

The second term represents the rotation of the 

supporting plate structure modeled here as an elastic 

spring. D is the flexural rigidity of the web plate, 
w 

From geometrical considerations: 
v = zB 

making a simple choice for B; 

B = B sin (m^ x/a) 
o 

substituting into the expression for the strain energy 
and integrating over a constant depth d results in: 
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V = D a d (rn 
w 



2 2 2 2 
/Va) B C d (m/7"/a) 
o 

2 

+ CB a/4 



o 



6(1 - V ) 3/12 



(3-1) 
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3. 2 Development of the Virtual. Work Eguat fon^. 

As for Model I, Model II is not assumed to be 

perfect but rather to have initial imperfect ions. These 

* * * 

initial imperfections v , w , B and the additional 
deformations v, w, and B are illustrated in figure 
(3-1). 




Figure (3-1) 1 1 lustrat ion of the initial imperfect ions 

v*, w*, B* and additional deformations v, w, B. 



For the case of an inplane axial load, the work 
equation is : 



W ■ 




From figure 



S u(y,z) 
x 

(3-1), it 



dy dz 

can be shown that : 
2 2 



- u ♦ 1/2 (v ) + l/2(w ) 



x 



X 



& a 
X 



X 
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and 



for a perfect web section. With imperfections, 

using the inextensibi 1 ity assumption with regard to the 

length of the stiffener setting & = 0, the result is : 

x 



*2 * 2 
u = 1/2 (v + v ) + l/2(w + w ) 

x xx xx 



u = 1/2 




* 2 *2 
C(v + v ) + l/2(w + w ) 3 dx 

xx xx 

(3-2) 



where u = u(y,z) 

As shown in figure (3-1), the geometrical 
relat ionships (equations 2-1) are used to simplify 
equation (3-2): 



f. 



a 2 2 



u (y, 2 ) = 1/2 (z + y ) (B + B ) dx 



A simple choice is made for B and B . 



B = B sin (m7hx/a) 
o 



* * 

B = B sin (m /7" x/a) 
o 



Thus: 

2 2 * 2 

u(y,z) = l/4(z + y ) (B + B MrnxVa) a 

o o 



The work equation now becomes: 



rr £ 2 * 2 
W = J / £ Cl/4 (z + y ) (B + B ) (m ^/a) al dy dz 

' J Ax o o 



46 



w = 



*2 2 

1/4 S (B + B ) <m 1? /a) a 
X o o 



I 

p 



3. 3 Determination of the Critical Buckling Stress.. 

Applying the calculus of variations with respect to 

B to the work and strain energy equations yields: 
o 

* 2 2 f 

d W = 1/2 S <B + B ) (m /?Va) al d B 

X o o p o 

r 2 2 2 
6 V = -C D ad (m^/a) B C d <m 7>/a) 



w 



+ 6(1 - V) 3/6 +CB a/2 > <f B 



Applying the principle of minimum potential energy: 

TF = S\? - Sw = 0 

<? W = & V 

Solving for % = £ : 

e CR 

2 2 2 

D d<m77Va> Cd (m^/a) + 6(1- V) 3/3 + C B 

w o 



CR 



(m » /a) 



* 

B + B 
o o 



(3-3) 

For the case of zero restraint against rotation; 
C = 0: 
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o 



2 2 2 

D d(m^ /a) Cd <m^/a) + 6(1-V)3 B 

w 

£ 

CR * 

3 1 B + B 

p o o 

(3-4) 

It is easily seen that the lowest £ for this case 

CR 

results from the simplest of mode shapes, rn=l. 

When the rotational restraint is not equal to zero; 

the lowest value of £ depends upon the degree of 

CR 

rotational restraint: 

3 1/4 

m = (a/7> ) (3C/D d ) (3-5) 

w 

With the use of equation (3-5), m may be evaluated 
arjd rounded up or down to the nearest integer for the 
correct wave shape. 

It is easily seen that as B becomes large £ 

o CR 

approaches a limit. Thus with (S ) = £ : 

CR LIMIT CL 



£ 

CL 



2 2 2 

D d (m ^/a) Cd (m7?Va> + 6(l-V)]/3 + C 

w 



2 

(m ^*/a) I 

P 



(3-6) 



Figure (3-2) illustrates the behavior of § with 

* CR 

increasing B + B . 

o o 

Further simpl if icat ion and rearranging of 
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. 



equations (3-3) and (3—6) provides: 



* 

B 

o 

B (3-7) 

o 

£ /£ - 1 
CL CR 




* 

Figure (3-2) Behavior of £ with B + B 

CR o o 



The first yield load will not be determined for 
this case under this method of analysis due to the 
relative complexity of equation (3-3) and due to the 
fact that the Perry-Robertson method will not lend 
itself to further simplification along this train of 
analysis. Since a simplified formula for design 
purposes is the goal of this thesis, the next chapter 
will use a different tactic to solve the problem. 
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4. MODEL II - GENERALIZED ANALYSIS 



4. 1 Method of Approach and Strain Energy Eguat fon^ 

The analysis thus far for the case of a flat bar 
stiffener under end loading has proceeded along the 
lines of thin plate theory up to the development of an 
expression for the critical buckling stress. 

This chapter uses beam theory to develop a simpler 
formulation in a manner analagous to that performed on 
Model I in chapter 2. In fact, this analysis will 
merely be a general i zat ion of the results for Model I 
applied to Model II. It must be remembered here that 
Model I consisted of two unequal flanges with the top 
flange deforming or rotating about the bottom flange. 
At no time in the derivation of the formula for the 
critical stress of the tee stiffener is there an 
assumption made about the position of the top flange 
relative to the bottom flange. Thus the same analysis 
with minor changes will be performed for the flat bar 
stiffener. However, in the case of the flat bar 
stiffener of Model II, there is a restraining effect 
provided by the juncture of the flat bar with the 
plate. This in effect can be modeled as a spring with a 
rotational constant C. Thus the strain energy equation 
would have four terms (instead of the three terms for 
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Model I). These terms represent (in order of 
appearance) sideways bending; longitudinal warping; 
torsion; and rotation of the supporting plate structure 
modeled as an elastic spring. Thus the strain energy 
equation is: 



V = 1/ 



2 X 



2 2 2 2 
[El v + EC B + GJB + CB 1 dx 

O Z XX W XX X 



Again choose: 

B = B sin(m 1^ x/a) 
o 

Substituting and integrating results in: 

2 2 4 2 

V = 1/4 a B C Ed s + C ) (m /a) + GJ (m Ti'/a) + C 1 

o Z W 

Applying the calculus of variations with respect to 



B yields: 
o 



<- 2 4 

SV = 1/2 a B C Ed s + C ) <m fiVa) 



w 



+ GJ (m T^/a) + C 3 cf 



B 



(4-1) 



4. 2 The Work Equation;. 

The assumptions used in deriving the work equation 
for Model I are the same as for Model II. Thus from 
chapter 2 we have: 



* 2 2 

W = 1/4 3 (B + B ) (m >7"/a) a I 
e o o p 



51 



and 



$ W = 



1/2 § (B 
e o 



B*> (m 7?Va> 2 
o 



a I 






B 



(4-2) 



4. 3 Det erminat ign of the Critical. Buckling Stress 



Using 


equations (4-1) and 


(4-2) and 


applying 


the 


principle 


of minimum potential 


energy and 


solving 


for 


the stress 


yields (with 5 = £ 


) : 







e CR 



2 4 2 

E(I s + C )(m7>*/a) + GJ<m7r/a> + C B 

z w o 

§ 

CR 2 * 

(m^/a) I B + B 

p o o 

(4-3) 

Let : 

2 4 2 

E(I s + C Xm7?Va> + GJ(m^/a) + C 
z w 

5 (4-4) 

CL 2 

( rn ^Va) I 

P 

Using (4-4) in equation (4-3) and solving for B 

< 

yields: 

* 

B 

o 

B = (4-5) 

o 

£ /S - 1 
CL CR 



Comparison of equation (4—3) with equation (2-1) 
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reveals the only difference is that equation (4-3) has 
an additional term in the numerator (specifically the 
spring constant C). 

Since it is very difficult to quantify the spring 
constant C and (hopefully) since its value is small in 
comparison to other terms in expression (4-3), C will 
be set to zero value and ignored for the remainder of 
this analysis. 

Thus for the case of C = 0, equation (4-3) reduces 
to equation (2-2). 



4. 4 Det ermination of the First Yi.el_d Load.. 

The determinat ion of first yield load for the 
generalized case proceeds in a similar fashion as for 
the case of Model I. 



Again 


the 


total 


torque consists 


of 


two parts 


(reference 


1>, 


St. 


Venant’s torque 


and 


a warping 



torque. However, in this case the warping torque is 
(reference 2) : 

T = -El d w 

w Z XXX 

and 

w = dB 
x x 

resulting in: 
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T = -El d B 

W Z XXX 

Thus the total torque is: 

2 

T = 6JB - El d B 

X Z XXX 

Choosing B and T: 

B = B sin (m 77 x/a) 
o 



T = T cos (m 7/ x/a) 
o 

Leaves: 

2 2 

T = ( ~P“ /a) CGJ + El d ( flVa) 3 B 
o z o 

Using equation (4-5) to eliminate B leaves: 

o 

* 

B 

2 2 o 

T = 07 /a) CGJ + El d ( 7T/ a ) 3 

o z 

(£ /£ 

CL CR 



(4-6) 



- 1 ) 



(4-7) 



Again sectional considerat ions can be used to 
evaluate the maximum stress. The maximum stress 
developed in the stiffener due to the torque is now: 

2 

(S > = -Ed B 

x MAXT xx 

Application of the above definition of B and a 

choice of x such that (£ ) is a maximum leaves: 

x MAXT 

2 2 

(S ) = Ed ( 7^/a) B 

x MAXT o 
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(4-8) 



Using equation (4-6) to eliminate B : 

o 

2 

Ed T ( V'/a) 

o 

(5 ) 

x MflXT 2 2 

CGJ + El d <Jt' /a) 3 



2 



Using equation 



(£ ) 
x MflXT 



(4-7) in equation (4-8) leads 

* 

B 

2 2 o 

Ed ( T /a) 



(§5 /S - 1) 
CL CR 



to: 



(4-9) 



flg ain, from sectional considerat ions, the total 
maximum stress can be described by: 



£ = (§ ) + S (4-10) 

y x MflXT CR 

Using equation (4-9) in equation (4-10), solving for 

2 

% and dividing by £ to place into nondimensional 
CR y 

form yields: 



2 2 2 * 




(4-11) 



Choosing the smallest root and simplifying the 
equation by the use of the reduced slenderness ratio 
(equation 2-15) but redefining parameter A as: 



55 



A 



(4-12) 



2 2 * 
( "Va) Ed B 

o 



% 

y 



Results in (equation 2-17): 



3 (1 + 1/L + A) 

CR 



N 



2 2 2 
(1 + 1/L + A) - 4/L 



£ 2 

y 

(4-13) 
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5. COMPARISONS WITH OTHER STUDIES. 



5. 1 Comgar ison between Thin Plate and Bearn Theory 
Results 



In chapter 3, the critical buckling stress of a 
flat bar stiffener was derived through the use of thin 
plate theory: 

2 2 2 

D d(m7?Va) Cd (rnTTVa) + 6(l-V)Il/3 + C B 
w o 



CR 



(m #/a) I 



* 

B + B 
o o 



„ (3-3) 

In chapter 4, the critical buckling stress of a 
flat bar stiffener was found using one dimensional beam 
theory: 

2 '4 2 

E(I s + C ) (m 7/Va) + GJ(m"Va> + C B 

2 W O 



5 



CR 



(m^/a) I 



# 

B + B 
o o 

(4-3) 
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It should be noted here that when: 



2 2 2 

D dtm/TVa) Cd (m 7?"/a) + 6 ( l - V )]/3 = 



w 

2 


4 2 



E(I s + C ) (rn Tr/ A ) + GJ(m7?Va) 
z w 

the agreement between the two methods would be exact. 

For the elastic range of experimentat ion and with: 
3 

Et 



w 

D = 

w 

12(1-^ ) 


G = E/2C1+V) = E/2.6 


3 

I = (d t )/12 
z w w 


3 

J = (d t >/3 
w w 


d = d 

w 


11 

Q 

• 

Ul 



The above relationship reduces to: 



2 2 

d (m TT /a) 1 

w 

+ ? 


2 2 

d (rnTT'/a) 1 

w 


SS 

10.92 2.6 


12 7. 8 


(Thin Plate Theory) 


(Beam Theory) 



From a cursory examination it would appear that the 



Thin Plate Theory results 


(equation 3-3) will probably 


give a larger value for 


the buckling stress for the 


normal range of stiffener 


dimensions used in ships. To 



test this hypothesis, the dimensions of a flat bar 
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stiffener used in experiment #21 of reference (6) will 

be used. 

where: 

d = 87. 4 cm 
w 

a = 730 cm 
m = 1 



For the thin plate theory: 



2 2 

d (mTr/a) 1 

w 

+ = 0. 398 

10.92 2.6 



For the beam theory: 

2 2 

d (m T/Va) 1 

w 

+ = 0. 14 



12 



7. a 



Thus, within the elastic range, the thin plate 
theory results would indeed provide higher buckling 
stress. The beam thoery results may thus be too 
conservative. However, section 5.3 shows good agreement 
between equation (4-11) (first yield load) and 
experimental yield results. 
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5. 2 Comoar isons with Published Formulae. 



Reference (3) performed a similar analysis to that 
performed in this thesis but without the inclusion of 
initial geometric imperfections. 

For the case of a Tee stiffener subjected to an 
axial load, equation 21 of reference (3) allows (in the 
notation of this thesis): 

2 2 2 
GJ + E(I s + C MmffVa) + C(a/mT?-) 
z w 

S = 

CR 

I 

P 

Equation (2-2) of this thesis, for a Tee stiffener 
in an analysis which ignores the contribution of the 
web and includes the contribut ions of initial defects, 
is: 

2 2 

[E(I 5 + C ) (m TTVa) + GJ3 B 



z w o 




I B + B 

p o o 



Thus equation (2-2), for the case of an initially 
perfect Tee stiffener reduces to equation 21 of 
reference (3) when the rotational restraint is held to 
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, 



it must be noted that since the 



be zero. However, 

physical models of the two derivations were slightly 

different - namely that the model used in this analysis 

consisted of two unequal flanges - some of the 

parameters of the equations will have different values 

( ie. I , J etc. > . 
z 

For the case of a flat bar stiffener, equation £2 
of reference <3) has: 



2 £ £ 
D dC(m^/a) d + 6(1-V)3 C(a/m/T) 

w 

£ = + 

CR 

31 I 

P P 



From chapter 4, equation (4-3): 

2 4 2 

E(I s + C ) <m 7/Va) + GJ(mr?Va) + C 

z w 

£ 

CR £ 

I (m 77-Va) 

P 



B 

o 



* 

B + B 
o o 



The following is a comparison for a Tee stiffener 
and a flat bar stiffener which was presented in 
reference (3). Figure (5-1) from reference (3) shows 
the geometric dimensions of the Tee and flat bar 
stiffeners. Table (5-1) shows how equations 21 and 22 
of reference (3) and equations (2-2) and (4-3) of this 
thesis compare with a finite element analysis performed 



61 



using a computer program developed and documented at 
the University of California (reference 12). 



4 

Flat Bar (S /E) x 10 C = 0 

CR 




FEM 


EGN (22) 


EGN (4-3) 


H 

II 

E 


6.44 


6. 42 


6. 23 


m = 2 


7.53 


7.51 


7. 03 


4 

Tee (5 /E) x 18 C = 0 

CR 




FEM 


EGN (21) 


EGN (2-2) 


m = 1 


26. 4 


26.3 


30. 69 


m = 2 


65.9 


72.6 


89.4 



Table (5—1) Comparison with Published Formulae 
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■ 0.635 cm - 



d * • 16 cm 



h ■ 2.42 cm 
A 



! 



t • 0 S4 cm 



b # • 30 cm 



FLAT BAR STIFFENER - • - 100 cm 




TEE STIFFENER - i - 160 cm 



Figure (5-1) Stiffener Geometries for 

Comparative Solutions 



S3 



The U. S. Navy has in the past used similar 
formula as those presented here for the deterrninat ion 
of critical buckling stresses. For the critical 
buckling stress of a Tee stiffener, in the notation of 
this thesis from reference 13 we have: 

2 

GJ + <mtt/a> EC 

w 

ej = 

CR 

I 

P 

where here: 

3 2 3 3 

C - f d t / 12 + d t 
w w f w 

3 2 3 

I =dt/3+dft +ft 
p w w w w f 



This is very similar to those results determined by 
reference (3), with the exception that the above 
formula does not account for the rotational restraint 
provided by the junction of the web and plate. Also I 

P 

is much simplified here and C , the torsion-bending 

w 2 

constant, equates to the quantity <1 s + C ) of 

2 w 

reference (3) - but also somewhat simplified. 
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5. 3 Comparisons with Exger imental R es ults. 



Reference (£) provides insight into the torsional 
buckling process through the experimental results 
published there. This report describes model tests on 
the collapse and post failure strength of Tee struts 
simulating flat bar stiffeners in a stiffened plate. 
The T struts were subjected to an axial load up to and 
beyond failure. The geometrical imperfections of the 
models were also recorded. The researchers of reference 
(6) (of Det norske Veritas) are to be lauded for the 
thoroughness and completeness of their report, 
particularly in the treatment of the experimental data. 

Table (5-2) contains the model data pertaining to 
this analysis for ten experiments of reference (£). 
Figures (5-2 to 5-11) show the load- deflection curves 
for these ten experiments. The solid line represents 
the experimental data. Point ft of each of these graphs 
represents the first yield load point as determined by 
equation (4-11). 

Notice that there is very good correlation between 
the experimental first yield load point and the results 
of equation (4-11). 
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# 

Model a h t B £ $ S 

wo y CR u 


21 


703 


87. 4 


4. 0 


1.67E-02 


256 


233 


230 


22 


703 


87. 3 


4. 0 


5. 6E-03 


271 


263 


246 


23 


703 


87.7 


4. 0 


1. 15E-02 


273 


257 


234 


24 


703 


87. 7 


4.0 


6. 0E-03 


271 


262 


240 


25 


703 


87.7 


4. 0 


8. CE-03 


272 


260 


230 


26 


703 


87. 7 


4. 0 


8. 0E-03 


256 


244 


243 


27 


703 


87. 4 


4. 0 


6. 0E-03 


266 


257 


272 


28 


703 


88. 0 


4. 0 


7. 99E-03 


256 


244 


241 


29 


703 


87.8 


4. 0 


4. 96E-03 


259 


252 


218 


30 


703 


88. 1 


4. 0 


1.26E-02 


271 


253 


230 



Not at ions 

a - length of model (cm) 
h — height of stiffener (cm) 

t - thickness of stiffener (cm) 

w 

* 

B — angle of rotation (radians) 
o 

2 

S - yield stress (N/mm ) 
y 

£ - first yield load stress (N/mm ) 

CR 

2 

£ - ultimate stress (experimental results) (N/mm ) 

u 



Table (5-2) Data of Model Tests 



l *2 




STROKE ( ) t 



MODEL 21 



Figure (5—2) Model Test #21 of Reference (6) 
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0.00 0.70 



l.<0 2.10 2. BO 3.50 4.20 4.00 S.6o o.iO 7.00 

STROKE in«) L. OO 



MODEL 7 2 



Figure (5-3) Model Test #22 of Reference (6) 
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t '1C 




MODEt 23 



Figure (5-4) Model Test #23 of Reference (6) 
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c os 




MODEL 24 



Figure (5-5) Model Test #24 of Reference (6) 



0W 
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r\ T) *j 



c e >2 




•0.00 0.70 1.40 2.10 2.80 3. SO 4.C0 4.90 S.bO o . 3-0 7.00 

STPOK£<nn> t on 



MODEL J5 



Figure (5-6) 



Model Test #25 of Reference (6) 



V 
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0.00 0.70 1.40 a. 10 2.80 3.50 4.20 4.90 5.60 6.30 7.00 

STROKE (llfl) £ on 

MODEL 28 



Figure (5-7) Model Test #26 of Reference (6) 
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c 




Figure <5-0> Model Test #27 of Reference <6> 
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* ►- il LJ 




STROKE! m) t 

MOOEl 28 



Figure (5—9) 



Model Test #28 of Reference (6) 
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i/, h- Q. u to cn 



E 02 
3.00 



ri 

/ 

n 

n 

i 

n 

n 



a.?v> 
2.40 
2.10 
1.80 
1 .50 
l .20 
0.90 
0 . £0 
0.30 
O.on 




0.00 0.70 1.40 2.10 2.80 3.50 4.20 4. DO S.oU £.30 7.0'J 

stroke <r.n> c 



MODEL 29 



Figure (5-10) Model Test #£9 of Reference (6) 



75 



v. »». m -s 



K C?c 




stroke < nn ) L 



MODEL 30 



Figure (5-11) Model Test #30 of Reference (6) 
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6. CONCLUSIONS 



This thesis utilized a design approach to tackle 
the problem of torsional buckling of stiffeners with a 
good deal of success. Integrated into all phases of 
this analysis was a considerat ion of the geometric 
imperfect ions of the stiffener - both Tee and flat bar. 
These structural members were analyzed as simplified 
sections and the Perry-Robertson approach used in other 
areas of engineering was successfully applied here. 

Chapter 5 showed that the derived results of this 
thesis reduced to currently published formulae for the 
case of the perfect stiffened. The contribution of this 
thesis is that the initial irnperf ect ions were taken 
into consideration, using a simple design approach, in 
the derivation of the results presented here. 

Also shown in chapter 5 was a generally good 
correlation between results of this thesis and axial 
load tests on flat bar stiffeners. Comparison of the 
Tee stiffener formula with experimental data, which 
measures and reports initial imperfect ions, needs to be 
performed. Also, comparisons should be made with other 
tests on flat bar stiffeners in more dimensional ranges 
to determine the range of validity for the results 
derived here. 
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Further work needs to be done in this area. 



Specifically, determinat ion of the critical buckling 
stress and first yield load due to lateral loadings. 
Also, a simplified method of handling the mode 
interactions of combined lateral and axial loadings is 
needed. In addition, it is an asset to understand from 
the work reported in references (6) and (15) that the 
falling path of the load-deflection curve (figure S-l) 
is represented by the tripping mechanism which is not 
covered in this analysis. 
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Figure (6-1) Experimental Load Deflection Curve 
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